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Abstract. We derive easy-to-handle approximations for 
the polar gap electric field due to inertial frame dragging 
as derived by Harding & Muslimov (1998). A formula valid 
for polar gap height comparable to the polar cap radius is 
presented. 
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1. Introduction 

Accurate models for E\\ — the electric field component 
parallel to the local magnetic field B above polar caps of 
rotating magnetized neutron stars — are of great theo- 
retical interest in the context of magnetospheric activity 
of pulsars. In the framework of space charge limited flow 
model (originally introduced by Arons & Scharlemann 
1979^ ), Harding & Muslimov ( gggg ) (hereafter HM98) de- 



rived -Ey including the general relativistic frame dragging 
effect, worked out by Muslimov & Tsygan ( 1992| ) (here- 
after MT92). HM98 considered the case with Eu screened 
at both a lower and an upper boundary of acceleration re- 
gion i.e. at the star surface and at a pair formation front. 
Since the full expression for En is too cumbersome for 
practical applications, HM98 offered simple analytic ex- 
pressions valid in various limiting cases. In this paper we 
revise and extend their results. In particular, we show that 
in the most frequently considered case, i.e. when the gap 
accelerator length is comparable to the size of the polar 
cap radius, the approximation derived by HM98 overesti- 
mates Eu by a ratio of the neutron star radius to the polar 
cap radius. 

2. Formulae for Eu 

In the following, h denotes the altitude above the neu- 
tron star surface and h c is the gap height above which a 



scaled with the half-opening angle of the polar magnetic 
flux tube 9(rj). 

A solution of Poisson's equation for the polar gap with 
h c -C R ns as derived by HM98 reads 



E 
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where, 

S i = £ fc/ffifcf) z <=,7i), 
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T(z, z c , 7) = - [ai {m - l)e 72 + a 2 (777 + l)e 1Z + 

+ ai(l - 7) - a 2 (l + 7)]/(oi + 02), (3) 

a a = ( 77?c + i) e -T««- 7 -i ) oa=7-l-(7»?c-l)e 7 * e (4) 

and 
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where ki and ki are the positive roots of the Bessel func- 
tions Jo and Ji, respectively, H (1)6(1) w 1, and x is a tilt 
angle between the rotation and the magnetic dipole axes. 
Values of Eq, 9q, e, and k are given by: 



En — B 



pc c 



pair formation front screens Eu (i.e. Eu 
Altitudes z = h/R nsl z c = h c /R ns and radial distances 
r] = 1 + z, rj c = 1 + z c scaled with the star radius R ns 
will also be used. The magnetic colatitude £ = 0/0(rj) is 



for /i > h c ). e=^^, 



MR'f ' 
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where B pc is the magnetic field strength at the magnetic 
pole, Q is the pulsar rotation frequency, and /, M are the 
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Fig. 1. The electric field E\\ in the limit h c <C R ns as a 
function of height h in units of r pc for several values of 
gap height h c = 0.03, 0.1, 0.3, 0.44, 1.0, 3.0, and 10.0 r pc . 
Thick solid lines present E^P of Eq. (|l|). Open circles mark 

approximation E? ] (Eq. §for the cases 0.03, 0.1, 0.3, 0.44, 

and 1.0 r pc ). The dot-dashed line is for i5|| 4 ^ (Eq. 13) and 



dashed lines are for E^ ] (Eq. [TJ, the cases h c /r pc = 0.3, 
0.44, and 1.0). The thin solid line indicates the case with 
no upper gap boundary (si 3) of Eq. ph. The vertical 



dashed line marks h = R ns . For h c = 0.44 
imations E^ and e[^ coincide. We assumed B„ 



pci 



approx- 
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G, P = 0.1 s, 
M = IAMq. 



0.2 rad, £, = 0.7, R n 
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moment of inertia, and mass of the neutron star, respec- 
tively. The function f{rf) reads 



f(r)) = -3s 3 [In (1 - x- 1 ) + (1 + (2x)- 1 )/x] 



(7) 



with x ~ T]/e. 

The electric field E^' as a function of height h is drawn 
in Fig. 1 for several values of gap height h c . It can be 
seen that En saturates in a twofold way: First, for a fixed 
h c it assumes a constant value at h exceeding the polar 
cap radius r pc (after a linear increase with h) . Second, for 
a fixed h -C h c it initially increases linearly with h c to 
become constant for h c ;> r pc . As noted by Harding & 
Muslimov this behaviour can be reproduced with simple 
approximations of Eq. (|TJ) for different regimes of validity. 

Since ki > 2.4 it follows form Eq. (H) that jiZ becomes 
smaller than 1 for h ^ r pc /3. Thus, in the limit h r pc /3 



and h c <C r pc /3 (or jiZ ^ I and jiZ c <C 1) the function 
T(z, z c , ji) may be approximated with 



F(z, z c , ji) ~ -7? 
which leads to 



1 I zz r 
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(HM98), where the relations 
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have been used (see eg. MT92). As can be seen in Fig 1, in 
practice, Eq. (||) (open circles) reproduces Eq. (Q) (thick 
solid line) for h c < r pc /3 (and h < h c ). 

Approximate behaviour of E\\ for h c ~ r pc can be 
determined by considering the opposite limiting case for 
the accelerator height: h c ;§> r pc /3. Assuming jiZ c 3> 1, 
h < (h c — r pc ), (and h c <C Rns) in Eq. (Ill) one obtains 



T{z, z c , 7i) ~ 7i (l - e 7iZ ) 



(11) 



which results in the same formula for E\\ as the one derived 
by MT92 for the case with no upper gap boundary (with 
the condition E\\ = fulfilled only at the star surface): 



J5„ ^-S^o^cosx^^-^tl 
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Nevertheless, for r pc <C /i c -C i? ns , Eq. ( [l2] ) reproduces 
(Eq. 0) almost over the entire acceleration length 
(see the thin solid line in Fig. 1) except its very final part, 
where h ^ (h c - r pc ). 

Taylor expansion of (Jill) reveals the linear increase of 



Eii with h for h S5 0.1r pc and the saturation above h 
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Fig. 2. Magnetic colatitude functions and /2(C) (see 
Eq. (|l5| ) for definition). 



pc 



^ h < (h c - r pc ) and 2r pc <, h c < i? ns . The 



for 

magnetic colatitude profiles 



, and / 2 (0=E? 



are shown in Fig. 2 and the relations 
8J (fciQ a ^l_6Ji(fe,0 



(15) 



- fc? Ji(fc,) 



£(l-0 (16) 



have been used in derivation of ([14j) (eg. MT92). Approx- 



imation (13) is shown in Fig. 1 as a dot-dashed line. 
For h c <C R ns the symmetry 



Ef\h)~E\, ll (h c -h) 



(1), 



(17) 



may easily be proven to hold. Therefore, to extend the 
validity of Eq. |[|) to h < h c it is useful to construct the 
approximation 



z 

1 

Zr 



(18) 



which works reasonably well for 0.5r pc ^ h c ^ r pc (dashed 
lines in Fig. 1). 

The expression derived by HM98 for the case h c ~ r pc 
clearly overestimates En by a factor 9q 1 w (rpc/i^ns)" 1 
( cf. eq. (A3) in HM98, also eq. (1) in Harding & Muslimov 
( |l998a|) , or a formula used by Zhang, Harding & Muslimov 
( pOOOj )). Moreover, unlike in eq. (A3) of HM98, E^ ] de- 
pends on the magnetic colatitude £ via /i(£)- O ur approx- 
imations are correct because they reproduce the "exact" 



formula for (Eq. @). Moreover, E^ ] of Eq. (||) is 

identical to the formula of MT92 for the case with no 



screening at h c (taken in the limit 7^2 <C 1). This identity 
should hold since for r pc < h c <C R nB both gap boundaries 
influence E\\ only within corresponding adjacent regions of 
height 



' pc • 



Since Eq. (|17|) holds for h c <C R ns , E\\ in the upper 
part of accelerator (r pc < h < h c ) may be approximated 
with 



Ef\h c 



-h) 



Ef(h c -h) 
for (h c 



PC 



<>h< hr 



E 



(5) 



(19) 



for r pc < h< (h c 



In the limit where h 
of Poisson's equation for the elongated polar gap reads 



3> r pc and h 3> r pc , a solution 



E 
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(MT92), where now 

fe ■ 

ll ~ e(r, c )r, c (l-e/7, c )i/2 ' 
~ fci 

'» ~ e(, )c ) J?c (i- £ /^)i/2 I 

.9(77) - 6{rj)5{rj)H{ri), 6{ v ) = Q {vfO)/f{v)) l/2 , and the 
functions S, and iJ arc defined in MT92 or HM98. Since 
7, :3> 1, the two summation terms which reflect the screen- 
ing effect of the pair formation front at h c contribute 
to En only at h very close to h c (see Fig. 3) so that 
li(Vc)(llc - V) ~ 1- Thus, for h c > r pc and r pc < ft < ft c 
one may simply use a formula derived in MT92 for the 
case with no screening at h c : 



E 



(9) 



—Ei 



3 K 



+ ^.g(ry)e(l-e 2 )sin X cos0L (22) 



For P 



0.1 s the approximation E^ 

proximation E { 
at h ~ 0.6r 



(9) 



matches the ap- 



' pc 



^ (which is valid for low altitudes (Eq. |12|) 
This altitude is much lower than esti- 



mated by Muslimov & Harding ( |1997[ ). Consequently, be- 
low ~ 0.6r r 



' pc cither E^ or E^ should be used. 
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Fig. 3. The run of .Em as a function of h/r pc in the limit 
h c ^> r pc . We choose h c = 2R ns (marked with dashed ver- 
tical line); the other parameters are assumed as in Fig. 1. 
Thick solid line presents E,| of Eq. J2(i| ) . It coincides with 

the case without the upper gap boundary (J5jj 9) of Eq. @), 
dotted) except for h ^> h c . The dot-dashed line presents 
E\p of Eq. (|l^). Thin solid line presents the approxima- 
tion Ef ] of Eq. (H). 



estimates of the height and width of the acceleration zone 
and of the maximum particle energy need to be revised. 

Since for an elongated polar gap (h c 3> r pc ) the upper 
gap boundary (pair formation front) influences the electric 
field only within the negligible (<C h c ) upper part of ac- 
celerator, the formulae derived by Muslimov and Tsygan 
(1992) for no upper boundary at h c may be used within 
the entire gap with no significant overestimates of electron 
energies. 
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3. Conclusions 

We find that for accelerator height h c approaching the 
polar cap radius r pc the accelerating electric field (for 
nonorthogonal rotators) may be approximated according 
to Eq. ©: 

E \\ = - 1 - 46 h ( X " Q h(0cos X Gauss, (23) 

where B\i = B pc /(10 12 G), P is the rotation period in 
seconds, h is in cm, for a star with M = 1AM®, R ns = 10 6 
cm. This result differs from the approximation derived in 
HM98 by a considerable period-dependent factor r pc /R ns . 
Moreover, unlike the HM98' formula it depends on the 
magnetic colatitude £ through the factor /i(£). It may 
be of particular importance for a possibility of limiting 
acceleration by the resonant inverse compton scattering 



(see Dyks & Rudak g000|) 



We emphasize that the corrections presented in this 
research note do not affect the results of the numerical 
calculations presented by Harding & Muslimov ( 1998] ) be- 



cause their numerical procedures include the exact expres- 
sions given by Eqs. ([!]) and (20). However, their analytic 



